Hydrodynamics at large baryon densities: Understanding proton vs. anti-proton V2 
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We study the importance of the initial state, baryon stopping and baryon number transport for the 
dynamical evolution of a strongly interacting system produced in heavy ion collisions. We employ a 
hybrid model, which combines the fluid dynamical evolution of the fireball with a transport treatment 
for the initial state and the final hadronic phase. We present results for collisions at beam energies 
from y/SNN = 7.7 to 200 GeV. We study various observables such as the centrality dependent freeze 
out parameters, the non-monotonic behavior of effective slope parameter parameter with particle 
mass as well as the apparent difference in particle and anti-particle elliptic flow. Our results are in 
reasonable agreement with the available data. We find that the propagation of the baryon-number 
current in the hydrodynamic evolution as well as the transport treatment of the hadronic phase are 
essential for reproducing the experimental data. 
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I. INTRODUCTION 

Over the last decades many experimental programs at 
the Brookhaven National Laboratory and CERN facili- 
ties have been devoted to finding signals of a new state 
of matter, the Quark Gluon Plasma (QGP), by means of 
relativistic heavy ion collisions. These experiments have 
produced a wealth of data including particle ratios and 
yields, transverse and longitudinal momentum spectra as 
well as the coefficients of a Fourier decomposition of the 
transverse flow patterns [H-tltij. Several observations at 
the Relativistic Heavy Ion Collider (RHIC) indicate that 
a strongly coupled QGP (sQGP) dominates the dynam- 
ical evolution [17H20| . In addition the energy dependen- 
cies of various observables, such as the K/ir ratio, already 
show anomalies at low SPS energies which might be re- 
lated to the onset of deconfinement and chiral symmetry 
restoration at lower energies [ill [l6[ . While these are in- 
triguing observations, experience has taught us that the 
interpretation of experimental results and their relation 
to the deconhnement phase transition is often ambiguous 
and extensive model studies are required to understand 
the numerous observables. 

In order to study the physics of heavy ion collisions 
many models have been developed which address specific 
aspects of the these reactions. For example, hadron 
ratios are well described by thermal models which 
employ a hadronic resonance gas at a fixed temperature 
and chemical potential for the description of particle 
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yields. The study of transverse and elliptic flow ob- 
servables require more complex models: at the lowest 
energies hadronic transport, including the effects of 
hadronic potentials/interactions ( see e.g. [H], [UHll]) 
are applied, while at the highest RHIC energies a 
description of the system in terms of fluid dynamics 
seems to be successful. Since it is desirable to obtain a 
more comprehensive picture of the whole dynamics of 
heavy ion reactions, various so called hybrid approaches 
have been developed during the last years [26|. In 
these models one commonly uses initial conditions that 
are calculated in a non equilibrium model which are 
followed by an ideal or viscous hydrodynamic evolution 
coupled to a Boltzmann equation for the freeze-out stage 
[27H35j . Alternative approaches couple a partonic phase 
to a hadronic transport model [36H38| . As each phase 
may contribute considerably to final state observables, 
the interpretation of experimental results and their 
relation to the deconfinement phase transition are often 
a difficult task, and the contributions of the different 
phases need to be evaluated thoroughly. 

The purpose of this paper is to study heavy ion col- 
lisions within the framework of one such hybrid model, 
the so called UrQMD hybrid model [39| with the aim to 
re-evaluate some of the commonly accepted interpreta- 
tions of various observables. Specifically we will address 
the centrality dependence of freeze out parameters, trans- 
verse momentum spectra as well as anti-particle elliptic 
flow at low energies. 

This paper is organized as follows: First we briefly in- 
troduce the hybrid model. Next we discuss particle pro- 
duction and transverse spectra obtained with the model. 
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In the final section we present our results for the elliptic 
flow of particles and anti-particles with an emphasis on 
the difference between proton and anti-proton flow. 

II. THE HYBRID MODEL 

The UrQMD hybrid model combines the advantages of 
a hadronic transport model with an intermediate hydro- 
dynamical stage for the hot and dense phase of a heavy 
ion collision. The UrQMD Model [HI, [22| (in its cascade 
mode) is used to calculate the initial state of a heavy 
ion collision for the hydrodynamical evolution 40] . This 
is done to account for the non-equilibrium dynamics in 
the very early stage of the collision. The coupling be- 
tween the UrQMD initial state and the hydrodynamical 
evolution happens at a time i s tart when the two Lorentz- 
contracted nuclei have passed through each other. 



To transfer all particles back into the UrQMD model, 
an approximate iso-eigentime transition is chosen (see 
f45( | for details). To this end we apply the Cooper Frye 
prescription [46| to individual transverse slices, of thick- 
ness Az = 0.2fm, at a time-like transition hypersurface. 
The transition time for a given slice is determined by 
the time when the energy density e in every cell of this 
slice has dropped below five times the nuclear ground 
state energy density, i.e. below ~ 730MeV/fm 3 . As a 
result we obtain a longitudinal iso-eigentime transition 
with an almost rapidity independent maximum switch- 
ing temperature for beam energies above ^/snn ~ 10 
GeV. This ensures that all cells have passed through a 
potential mixed phase of the equation of state and the ef- 
fective degrees of freedom at the transition are hadronic. 
In a given slice the hydrodynamic fields are transformed 
to particle degrees of freedom via the Cooper-Frye equa- 
tion: 



tstart — 1 ry — W) 

where 7c. m. is the center of mass frame Lorentz fac- 
tor and R is the radius of the nucleus. At this start- 
ing time all initial collisions have happened. It is further 
the earliest time at which local thermodynamical equilib- 
rium may be achieved. At this time the energy, baryon 
number and momenta of all particles within UrQMD arc 
mapped onto the spatial grid of the hydrodynamic model 
by representing each hadron by a Gaussian of finite width 
(7 = 1 fm. In this approach the effects of event-by-event 
fluctuations and stopping of energy and baryon number 
density in the initial state are naturally included. 

The full (3+1) dimensional ideal hydrodynamic evolu- 
tion is performed using the SHASTA algorithm [U |42j ■ 
We solve the equations for the conservation of energy and 
momentum: 



on an isochronous time-like hypersurface in the 
computational frame (the hypersurface normal is 
da^ = (d 3 x, 0, 0, 0)). Here fi(x,p) are the distribution 
functions and §i the degeneracy factors for the different 
particle species i, which in our case are given by the 
appropriate Bose- or Fermi-distributions. 

The above transition procedure conserves baryon num- 
ber, electric charge, the total net strangeness and the 
total energy on an event- by- event basis. After the parti- 
cles are created they evolve according to a hadronic cas- 
cade (UrQMD) where final re-scatterings and decays are 
calculated until all interactions cease and the system de- 
couples. A more detailed description of the hybrid model 
including parameter tests and results can be found in [3t| . 



d^ u = (2) 

and for the conservation of the baryonic current N^, 

d^N" = (3) 

Here T^ u is the relativistic energy momentum tensor. 
The above partial differential equations are solved on a 
three-dimensional spatial Eulerian grid with fixed posi- 
tion and size 5x = 0.2 fm in the computational frame. 
To close the set of equations an equation of state (EoS), 
the pressure as function of energy and baryon number 
density p(e,n) needs to be specified. In the following we 
will use an EoS that corresponds to an hadron resonance 
gas. This is convenient as it includes the same degrees of 
freedom as the UrQMD model, which essentially allows 
us to study the effect of local equilibrium on the differ- 
ent observables. Previous investigations have shown that 
different equations of state only lead to insignificant dif- 
ferences in the results [43|, HiJ . 



A. Particle yields and the Temperature 
distribution 

The description of particle production in heavy ion col- 
lisions using statistical methods, such as thermal models 
47 49] and fluid dynamical models with a sudden freeze 
out, are usually based on the assumption that for a given 
incident beam energy chemical freeze out occurs at a fixed 
temperature and chemical potential. In a fluid dynamical 
treatment of the expansion this implies that for a given 
freeze out criterion the particle abundances are fixed by 
the total entropy per baryon (S/A) produced in the very 
early stage of the collision, since the subsequent expan- 
sion is isentropic. 

In all models used to describe the properties of the 
early stage of an heavy ion collision (Glauber model, ge- 
ometrical overlap, hadronic transport) the initial state is 
defined by the total energy and baryon number deposited 
in the fireball. In a Glauber model [l], [l(J H3] the energy 
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deposition in the transverse plane is proportional to the 
number of wounded nuclcons: 



e(x,y,T ) = 
KlT A (x+k,y) 



(•5) 
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where Ta and Tb are the nuclear thickness functions of 
the incoming nuclei A and B, a is an energy dependent 
cross section, and if is a proportionality factor. If we 
further assume that the initial baryon density is propor- 
tional to the initial energy density (To| , 



n{x,y,T ) = C(y?s NN ) ■ e(x,y,T ) , 



(6) 



the initial energy per baryon (E/A)(r) as a function of 
the transverse radius is a constant which only depends 
on the collision energy. At the same time the energy 
and baryon densities drop rather quickly as a function 
of radius r. Since E/A is constant, a decreasing energy 
density implies that the entropy per baryon S/A must 
vary as a function of r as the entropy density in general 
does not scale with the energy density, 



s(x,y;T ) 
e(x,y;r ) 



constant. 



(7) 



In addition the average energy density varies as a func- 
tion of centrality. Since E/A is constant we also expect 
a centrality dependence of the total entropy per baryon. 

Figure Q] exemplifies the effects of the initial geometry 
on the final chemical composition of the fireball. The 
upper plot of figure Q] shows the total entropy per 
baryon produced in collisions of Au nuclei at energies 
of \J (snn) = 7.7-19.6 GeV, as a function of impact 
parameter b. The entropy per baryon is calculated, av- 
eraging over 1000 UrQMD initial state events (although 
it has been shown that S/A does only vary weakly on 
an event-by-event basis [40]). A Hadron Resonance 
Gas (HRG) equation of state is used to estimate the 
values for the entropy per baryon as it best resembles 
the active degrees of freedom at the particle freeze 
out. We find that the produced S/A increases with 
the impact parameter, simply because of the different 
initial geometry. In the lower part of figure [1] we also 
present estimates for the influence of this effect on 
the baryon chemical potential fiB- Using the HRG 
we can determine the baryonic chemical potential /j,b, 
corresponding the value of S/A at a fixed temperature 
Tf, as a function of b. It is clear that fiB decreases 
with increasing b. We note that such an effect has 
been observed in experiment (5l| . where thermal fits 
to different centrality selections show a decreasing /j,b 
with centrality which is of comparable magnitude to our 
results shown here. This observation indicates that the 
systems created in heavy ion collisions of fixed energy 
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FIG. 1: (Color online) Upper plot: Centrality dependence of 
the Entropy per baryon produced. Lower plot: Corresponding 
baryon chemical potential at fixed temperatures Tj. 



but varying centrality cannot be characterized by single 
values for the thermal parameters. 

An important aspect of particle production in hybrid 
models is the treatment of the transition from the hydro- 
dynamical phase to the the transport phase. As we have 
discussed above (see also we transform full trans- 
verse slices at a given time i, using the Cooper Frye pre- 
scription. This implies that contrary to an isothermal 
transition particles are emitted over a range of temper- 
atures/densities. While this may no be a standard pro- 
cedure it has been shown in several publications (39l. |52|- 
|56| that this hybrid model is very well suited to describe 
particle ratios and yields over a wide range of energies. 
In addition, recent comparisons of an iso-thermal 3+1 di- 
mensional hypersurface [57J indicate that deviations only 
occur at large transverse momenta (p± > 1 GeV), corre- 
sponding in our freeze out prescription to very cold and 
fast cells which are being treated as equilibrated for an 
extended period of time. Since all observables we discuss 
in this paper are mainly determined at lower momenta 
our results hold also for the iso-thermal case. The com- 
parison also implies that both dynamical descriptions, 
transport and hydrodynamics, give equivalent results for 
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FIG. 2: (Color online) Probability distributions P(T) for pi- 
ons (black solid line) and O - 's (red dashed line) to be pro- 
duced from a cell with a give temperature T. The results 
shown are for central collisions of Pb nuclei with a beam en- 
ergy of y/SNN = 17.3 GeV. 



the fireballs expansion in the transition region consid- 
ered. 

Figure [2] shows the resulting probability distributions 
P(T) for pions and Q~'s to be produced from a cell 
with a given temperature T. The results shown are for 
central collisions of Pb nuclei at a center of mass en- 
ergy of \J sjv./v = 17.3 GeV. One can clearly see that 
the pions, due to their smaller mass, are generally emit- 
ted/produced at lower temperatures than the Si's. In 
other words, because we allow for emission from all lo- 
cal temperatures the heavier particles are preferentially 
emitted from regions of higher temperatures. The lower 
plot of Fig. [3] illustrates how the average value of Tcf 
changes as a function of the particle mass for different 
hadron species. Results for Pb+Pb/Au+Au collisions at 
^snn = 17.3 and 200 GeV are shown. One observes a 
clear trend of increasing Tcf with the particle mass, as 
shown explicitly for the pion and fi. 

We also observe a mass dependence of the average 
transverse flow (v±), defined as the average hydrody- 
namical radial flow velocity at the production points of 
the particles. Results for the different particle species 
are depicted in the upper part of figure |3] It is clear that 
hot and dense cells usually have a smaller flow velocity 
because they tend to be located at the center of the 
collision zone. The colder cells are more peripheral and 
therefore have acquired more flow. 

At lower beam energies, the finite baryon density 
(chemical potential) plays an important role. Because 
the proton has a large chemical potential, it is much 
more abundant at lower temperatures than it's anti- 
particle, resulting in a distinctly different value of (Tcf}- 
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FIG. 3: (Color online) The average transverse flow velocity 
(upper figure) and Tcf (lower figure) at the Cooper 
Frye transition, as a function of the particle mass for differ- 
ent hadron species. Results for most central Pb+Pb/Au+Au 
collisions at ^/snn = 17.3 (black squares) and 200 GeV (red 
circles) are shown. Note the anti-proton mass has been shifted 
for visibility. 



Figure |4] depicts the averaged ratio of anti-protons 
over protons, taken at the last time step of an Au+Au 
collision at a beam energy of ^snn — 11 -5 GeV, 
p/p = exp (— 2 (ib/T), as a function of radius (red 
dashed line). Clearly the p/p ratio is largest in the 
center of the collision zone where the radial flow has 
its minimum. The transverse flow velocity v± (black 
solid line) also extracted from our hydrodynamical 
calculation, increases linearly with the distance from the 
center of the collision. Consequently the anti-protons 
will acquire a smaller average transverse flow as com- 
pared to the protons. This conclusion also holds true 
for an iso-thermal transition because the ratio p/p will 
in general not be constant over the hypersurface. Note 
that the same argument should also hold for other 
conserved charges, such as strangeness and the third 
component of the isospin. Because the colliding nuclei 
have more neutrons than protons the third component 
of the isospin is finite, and negative in the produced 
fireball. Therefore, we expect the ratio tt + /tt^ to 
depend on the radius similar to the p/p-ratio, which 
in consequence will lead to a different transverse flow 
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FIG. 4: (Color online) Average transverse hydrodynamical 
flow velocity (black solid line) and the scaled ratio of pro- 
tons over anti-protons w exp (— 2/is /T) (red dashed line) as 
a function of transverse radius. The results shown are for 
the central transverse plane at the end of a hydrodynamical 
calculation for Au+Au at y^ijvjv = 11.5 GeV. 



velocity of ir + compared to 7r~. However the effect 
will be much weaker for pions since it depends on the 
ratio of the fugacities exp (— 2/j.i/T) and the isospin 
chemical potential is considerably smaller then the 
baryon chemical potential < /is < ^b- For the top 
SPS energy a thermal model analysis found values of: 
uj_= -5.0 MeV, us = 71.1 MeV and \i B = 266 ± 5 MeV 
581. 



III. THE SLOPE PARAMETER 

After discussing the integrated yields, let us next dis- 
cuss the transverse momentum spectra, which are consid- 
ered to be more sensitive to the decoupling temperature, 
Tth, and transverse flow than average yields. Taking Eq. 
(j4]) and neglecting any flow the transverse momentum 
distribution (for p z —¥ 0) is given by 



1 



dNi 



mj_ dm±dy 



cxm i (exp((m ± -/i l )/T)±l)- 1 (8) 



with the transverse mass mj_ = yj (m? + p\). The neg- 
ative slope of log (^3- d ^ dy ) gives an effective tempera- 
ture T e jf 63, 64], which, however, does not correspond to 
the actual temperature since the spectra are blue shifted 
due to transverse flow. To account for the flow the effec- 
tive temperature can be related to the decoupling tem- 
perature by T eff = T t h + mi (v±) 2 . 

These type of fits to experimental data typically result 
in a non-monotonic behavior in T e ft as a function of 
particle mass, which cannot be accounted for with a fixed 
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FIG. 5: (Color online) Transverse mass distributions of pions 
protons and fi's at mid rapidity \y\ < 0.5, divided by m\, for 
most central collisions of Pb nuclei at ^/sjvjv = 17.3 GeV. We 
compare the distributions from the hybrid model extracted 
at the final state (black squares) with those directly at the 
Cooper Frye transition with (red circles) and without (green 
triangles) resonance decays. Data from the NA49 experiment 
are depicted as blue circles [60l462| . 



temperature freeze out [65J, |66| . 

In addition to flow effects resonance decays as well as 
scattering processes in the hadronic phase affect the ex- 
tracted value of T e ff as a function of mass [H, HH, [63, [fi8[ . 
To successfully interpret the experimental results one, 
therefore, has to disentangle the different contributions 
to the transverse spectra in a consistent approach. 

Using the UrQMD hybrid model we can disentangle all 
the important contributions and explore what informa- 
tion about the hot and dense phase can be extracted from 
the spectra. Figure [5] shows the mid rapidity (\y\ < 0.5) 
m± spectra of pions protons and fi's, divided by m^_, for 
the most central collisions of Pb nuclei at \/(snn) = 17.3 
GeV. We compare the final spectra (black squares) with 
the ones obtained directly after the Cooper Frye transi- 
tion where we either let all resonances decay at the tran- 
sition (red circles) or do not allow for resonance decays at 
the transition, (green triangles). For comparison we also 
show the experimental data by the NA49 collaboration 
[6oTj62| as blue circles. 

We see that the spectra of both pions and protons are 
significantly modified due to hadronic interactions and 
resonance decays. We further find deviations of the spec- 
tra from an exponential shape already at the Cooper Frye 
transition, when no resonance decays rare taken into ac- 
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FIG. 6: (Color online) Ratios of transverse mass distribu- 
tions after to before resonance decays, shown for pions pro- 
tons kaons and fi's, as a function of m±. 



count. This is simply a result of the blue shift due to 
the finite flow already present at the transition, as shown 
in Fig. [3J We note that similar modifications of the p± 
spectra due to final state interactions were obtained in 
an earlier calculation [65[ , where the UrQMD model was 
used for the final state of a hydrodynamical calculation 
with an isothermal freeze out. 

In both calculations, ours as well as that of Ref. .655, 
only the fi does not change noticeably during the after- 
burner phase because of it's small hadronic cross section. 
Furthermore there are no resonance states in the model 
that contribute to the f2 spectrum. 

Resonances change the momentum spectra mainly 
at low momenta, because of the restricted phase space 
of the decays and their effects should become negli- 
gible at some point. This is demonstrated in Figure 
IHl where we show the ratio of the transverse mass 
spectra after the resonance decays to those before 
the decays for four different particles. The strongest 
modification of the spectra is observed for the pions 
and protons: they are enhanced by a factor of three 
at low momenta. However, we find that the effect of 
resonances is still considerable even for m± — mn > 1 
GeV, which is consistent with previous works [fM [67| . 
Therefore, it is questionable that an exponential fit 
below a transverse mass of mi — mo < 1 GeV is in any 
way justified, and even at higher transverse masses one 
would expect small deviations from the exponential form. 

In Fig. [7] we show the exponential fits to spectra 



FIG. 7: (Color online) Extracted slope parameters T e // from 
the hybrid model for different particle species. We compare 
final state results (full symbols) with results obtained at the 
cooper Frye transition without resonance decays (open sym- 
bols). 



obtained with the hybrid model over a mass range of 
I < m_i_ — m < 2 GeV. The filled symbols correspond to 
fits of final state spectra, including resonance decays and 
the UrQMD afterburner, for most central PbPb/AuAu 
collisions at \/(snn) = 17.3 GeV (black squares) and 
\J {snn) = 200 GeV (green diamonds) GeV. The red 
open circles correspond to fits to spectra directly after the 
Cooper Frye transition without resonance decays. One 
clearly observes a non-monotonic behavior in the effec- 
tive slope parameter as a function of the particles mass. 
It was shown in earlier publications [65[ that such a be- 
havior can also be accounted for when UrQMD is used 
as the afterburner for an iso-thermal transition. In that 
work the scattering processes in the transport phase cre- 
ate a situation where the actual dynamical decoupling 
takes place over a range of densities/temperatures lead- 
ing to results comparable with our study. As indicated 
before, the hydrodynamical and transport description 
give equivalent results over a certain range of densities 
and our non-isothermal transition is generally not ruled 
out. The question whether there is a sharp transition, as 
in the iso-thermal case, or a certain transition region may 
be further investigated by carefully studying the spectra 
of the </> and f2, who do not re-scatter much in the trans- 
port phase. 

To estimate how well the approximation 
T e ff = T t h + rrii(v±) 2 , with (v±) taken from fig- 
ure [31 compares with our fit results, we plotted the 
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simple estimates of T e ft as black crosses. They should 
be compared with the open red circles. It is obvious that 
the simple formula fails to describe the slope parameter 
for most of the light mesons and works more reliably for 
the heavy baryons. 



created in a heavy ion collisions dHl H, [z3| El llO] ■ 
The elliptic flow parameter v 2 is defined as the coefficient 
of the second Fourier component of the azimuthal distri- 
bution of the emitted particles: 



IV. ELLIPTIC FLOW 

The elliptic flow parameter v 2 has long been proposed 
as a probe for the properties of the hot and dense system 

I 



dN 
p±dp±d$> 



^- m [1 + 2v 2 {pjl) cos(2$) + 2v i (p±) cos(4$) - 

Z7T p±dp± 

I 



(9) 



where the azimuthal angle $ = (j)—ipRp is measured with 
respect to the direction of the reaction plane, ipRp 

The average or integrated elliptic flow coefficient, v 2 , 
is defined as: 

v 2 = (coa[2($)]) (10) 
and is given in terms of the azimuthal distribution as 

t , r±dS cos(2<E>) 



Jd P± J "d$ T 



'•2 = — — ^— — (n) 

(12) 
(13) 

The integrated elliptic flow is related to the p± dependent 
elliptic flow by 




V 2 = 



Jd P± v 2 ( P± )^L 



(14) 



To calculate v 2 in fluid dynamics one usually performs 
a Cooper-Frye freeze out and then uses the above defini- 
tions. 

Alternatively one can use the energy momentum tensor 

T, 

hydrodynamical computation [72j, |7- 



^ u and directly extract the flow anisotropy from the 



«2,T 



y^xx -^3 



mil 



f d 3 r f d 3 p &J&f {r>p) 



Jd*rfcPp&g&f(r,p) 



(15) 
(16) 



Here, f(r,p) is the phase space distribution, and 
T xx , T yy are the appropriate components of the energy 



momentum tensor. This prescription is different from 
the actual definition of v 2 which is given by Eq.(fT2"j). 
It has an additional weighting of p]_/po in front of 
dpfd® ■ However, it was empirically found [72j that the 
pion elliptic flow can be related to v 2 ^t as 2 v 2 v 2t T- 
Assuming that v 2 (p±_) oc p± and that f(r,p) is repre- 
sented by a Boltzmann distribution on can easily show 
that | v 2 ~ v 2 .t- However, the distribution function 
for pions is considerably distorted by resonance decays 
and far from being of Boltzmann type, especially at 
small transverse momenta. This brings the factor re- 
lating the two definitions close to the observed value of 2. 



In figure [S] we show results for the integrated v 2 of dif- 
ferent particle species, pions protons and anti-protons, 
as a function of time, extracted from the hydrodynam- 
ical phase of the hybrid model. For this calculation we 
used a non-fluctuating initial condition created by aver- 
aging 1000 UrQMD events of collisions of Au nuclei at an 
\J snn — 11-5 GeV and an impact parameter of b = 8fm. 
The symbols denote values of v 2 extracted from sampling 
the Cooper-Frye equation on an isochronous hypersurface 
at each time step, hence representing the correct defini- 
tion of v 2 . The lines represent values of v 2 .t extracted 
from the hydrodynamical energy momentum tensor as 
described in equation (|15p . For the pion elliptic flow we 
used the full energy momentum tensor and multiplied 
v 2 .t by a factor of 0.5 as suggested in [72j and one can 
observe a very good agreement of the methods. To ex- 
tract the proton and anti-proton flow we used the partial 
energy momentum tensor of the protons and anti-protons 
and multiplied them by 2/3 as suggested above. The par- 
tial 's for protons and anti- protons are related simply 
by: 



exp(-2fi B (x,y)/T(x,y)) 



rpp 



(17) 
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FIG. 8: (Color online) Elliptic flow as a function of time for an 
averaged UrQMD initial condition of two Au nuclei colliding 
at y^Sjvjv) = 11. 5 A GeV and an impact parameter of b — 8 
fm. The lines correspond to V2 extracted using the energy 
momentum tensor Eq |15l and the symbols denote results from 
a sampling of the Cooper Frye equation. 



FIG. 9: (Color online) Difference of particle and anti particle 
V2 as a function of beam energy for different particle species. 
The open symbols denote results at the Cooper-Frye tran- 
sition while the full symbols represent results with the full 
UrQMD afterburner stage. 



(pB and T being the baryon chemical potential and 
temperature respectively) . 

The resulting proton and anti-proton iVs are in reason- 
able qualitative agreement for both methods presented. 

In either case we observe that the integrated elliptic 
flow of protons is systematically larger than that of the 
anti-protons. This difference is also observed in recent 
preliminary experimental data from the STAR Collabo- 
ration [74|]. Below a collision energy of ^snn < 60 GeV 
the measured elliptic flow of particles is considerably dif- 
ferent from that of their anti-particles and the difference 
increases with decreasing beam energy and correspond- 
ingly increasing net baryon density. A recent transport 
calculation was able to explain such an effect qualita- 
tively by the inclusion of mean field type nuclear poten- 
tials [73]. Our calculations, however, indicate that the 
average flow of e.g. anti-protons is different from that 
of protons simply because their local 'weight', given by 
exp (— 2(1b(%, y)/T(x, y)) in the evaluation of e.g. equa- 
tion ([17} , varies due to the finite net baryon density. As 
shown in Fig. 2] the ratio of p/p is large for the cells in 
the center of the collision which have a high temperature 
and small flow viscosity, and it is small for the colder 
cell at the surface, which carry high flow velocities. As 
a result the average transverse flow for protons is larger 
than that for anti-protons resulting in a larger value for 



V2 of protons. 

Let us next quantify the discussed difference in ellip- 
tic flow and investigate to which extent our result is 
modified when we apply the full hybrid model includ- 
ing event-by-event fluctuations, resonance decays as well 
as an afterburner stage. Figure [9] shows the difference of 
particle— anti-particle V2 of different particle species as a 
function of ^snn for collisions of Au+Au nuclei at an 
impact parameter of b = 8 fm. In the plot the full sym- 
bols correspond to results obtained after the UrQMD fi- 
nal state while the open symbols denote results after the 
Cooper-Frye transition, including resonance decays. In 
general the difference in elliptic flow is most pronounced 
directly after the Cooper Frye transition out of the hy- 
drodynamical evolution and is washed out by the subse- 
quent UrQMD transport phase. For baryons we observe 
a considerable increase in the difference of V2 between 
particles and anti-particles with decreasing beam energy. 
The value of V2 is essentially the same for tt + and 7r~ . For 
the Kaons we observe a trend, which is opposite to the 
one found in data [74j ]. showing that K~ V2 seems larger 
than that of K + . However, this is only the case after the 
final UrQMD transport stage. Directly after the hydro- 
dynamical evolution kaons appear to have the same flow. 
The observed phenomenon can therefore be interpreted 
as a natural result of a net-density and a chemically equi- 
librated phase. 
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FIG. 10: (Color online) Effect of the final hadronic state on 
particle and anti-particle V2 as a function of beam energy. A«2 
is denned as the gain of the elliptic flow coefficient obtained 
in the final transport phase: vp nal - yCoaper-Frye_ 



FIG. 11: (Color online) Elliptic flow of protons, anti-protons 
and pions as a function of p±. We show results from an av- 
eraged UrQMD initial state (lines) as well as from the full 
hybrid model calculations after the Cooper Frye transition 
(No final state interaction) (symbols). 



Many transport descriptions fail to even qualitatively 
describe the phenomenon because they lack certain 
interaction channels which are important at the en- 
ergies considered here. For example the difference in 
the Kaon elliptic flow could be explained by missing 
strangeness exchange reactions in the final state. For the 
anti-particles the inclusion of all pair creation processes 
is important. Since anti-particles are very rare in low 
energy collisions, their flow is sensitive to their explicit 
interactions, i.e. annihilation and recreation, in the 
transport phase. At the highest beam energies a large 
number of anti-particles is produced at hadronization 
and regeneration has only a small effect on their 
abundance (seen in LHC data (7f| [77||). However when 
the beam energy is decreased, only few anti-particles 
are produced and the pair creation due to multi pion 
reactions becomes non negligible for the anti-particles 
bulk properties [78148 1|. 

To illustrate the effect of the hadronic transport 
phase on the elliptic flow, in Figure [TU] we show the 
difference between the value of vi in the final state to 
that obtained right after the Copper Frye transition 
for protons, anti-protons and charged pions. Clearly, 
at low energies the value for V2 for anti-protons in- 
creases appreciably during the hadronic transport phase. 
This is simply an effect of the annihilation process. 
Anti-protons moving in the out-of-plane (y) direction 
encounter more protons to annihilate with than those 



moving in the in-plane (x) direction. Since there are 
many more protons than anti-protons, the annihilation 
only affects the anti-protons. Because the reverse 
process nir — > p + p is not included in the transport 
model it is not clear how meaningful the final state 
effects (in the UrQMD phase) are for the anti-particle v-x- 

Note that in our calculation the difference between 
proton and anti-proton elliptic flow is caused solely by 
the non-zero net baryon number density, and chemical 
potential, and therefore a similar effects should be 
be observed for a finite isospin and net strangeness 
chemical potential at mid rapidity [82j if included 
properly in the model. As discussed at the end of 
section |TT] the difference in the (elliptic) flow of particles 
and anti-particles is due to the different local weighting 
cxp (— 2/itj(x, y)/T(x, y)) when evaluated over the hyper- 
surface. We can assume that the chemical potentials 
Hi have a different magnitude —fJLi < fJ,s < Ms an d 
therefore we expect the effect becomes weaker for Kaons 
and even more for pions. This is in fact the case for the 
experimental data, as tt~ show more vi than 7r + and K + 
more than K~, where the difference is smallest for the 
pions. Since the isospin chemical potential is negative 
due to the higher abundance of neutron in the colliding 
nuclei, this also explains why, contrary to the protons, 
in case of the pions the elliptic flow of the anti-particle, 
i.e. 7T~, is larger than that of the particle (tt + ). 
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FIG. 12: (Color online) Integrated elliptic flow of protons and 
anti-protons for any given definition of the Event Plane angle 



In figure [3] we see that, during the hydrodynamical 
evolution, the average transverse momentum of p is 
smaller than that of p. Therefore, part of the difference 
in the integrated r>i can be understood as weighting 
^2(p_l) at higher values of p± in case of the protons 
(see Eq. (fT4|) ). On the other hand the preliminary data 
by STAR indicate that also the differential elliptic 
flow parameter, t>2(pj_), differs between protons and 
anti-protons. To investigate this in Fig. [TT] we compare 
the differential elliptic flow for protons, anti-protons 
and pions after the Cooper-Frye transition. The lines 
represent results with the averaged initial conditions 
as described above, where we sampled the freeze out 
hypersurface to obtain sufficient statistics. We indeed 
observe a small difference for V2{pi_) between p and p . 
From the event-by-event calculations we obtain results 
depicted as symbols. Here we also only observe a small 
difference, within the still considerable errors. In any 
case, in our calculation we find the effect due to the 
different (p±) to be dominant. 

Another effect which can artificially increase the ob- 
served asymmetry of baryon vs. anti-baryon elliptic flow, 
and in particular the difference in p± dependent V2 , may 
arise from the way i>2 is measured in experiment. While 
in our calculation the reaction plane is well defined, ex- 
periments have to reconstruct the so called Event Plane 
in order to infer the Reaction Plane angle ^pp. This 
analysis method involves correlating the azimuthal angle 



4> of each particle with an event plane angle $ep which 
is defined by the other particles of that event 83]. Due 
to fluctuations, the Event Plane is usually not perfectly 
aligned with the Reaction Plane (see e.g. [Hj]). 

Furthermore, as a result of baryon number fluctuations 
an event plane direction, ^^ P , defined only by protons 
may deviate from that defined solely by anti-protons, 
\l/^p. As a consequence, the elliptic flow for anti- protons 
calculated with respect to the event plane of the protons 
would be slightly smaller than if it were calculated with 
respect to the event plane of the anti-protons. To illus- 
trate this we have plotted in Fig. [TT] the integrated V2 for 
protons and anti-protons as a function of a trial event 
plane angle 'FepJUI]. Obviously, for ^>ep = ^ep an< ^ 

ep — ^ep the values of V2 are maximal for protons 
and anti-protons, respectively, indicating the correct di- 
rection of the respective event planes. Furthermore, one 
can clearly see that the maxima for protons and anti- 
protons are separated by a finite angle. 

Since at low collision energies protons are much more 
abundant than anti-protons, it is likely that the determi- 
nation of the event plane is biased towards the direction 
of the proton Event Plane, 'f'fp. As a consequence, in 
each event the value extracted for V2 of anti-protons will 
be less than the maximum, as shown in Fig. 111! Thus, 
the event averaged value for the integrated elliptic flow 
for anti-protons will be systematically smaller than that 
of protons, even if their true values would be the same. 

At (nearly) vanishing net baryon densities this is not 
of importance because one has as many baryons as anti- 
baryons and there is no bias towards either the proton 
or anti-proton event plane. In the case of baryon anti- 
baryon asymmetry however, the Event Plane definition 
will always be biased towards the particle plane. In the 
energy range considered we estimated this effect to con- 
tribute to the difference of the flow by about 5 — 10 %, 
by averaging the relative difference of the peak positions 
over many events. Note however, that this effect could be 
excluded if the Event Plane/Reaction Plane can be mea- 
sured independently, e.g. from the spectator fragments, 
thus eliminating any bias. 



V. SUMMARY 

We have presented results calculated within the 
URQMD hybrid model for various observables from 
heavy ion collisions. The hybrid model is able to take 
into account the initial stopping, the explicit propagation 
of the baryon number, and the non-equilibrium transport 
in the hadronic phase, all of which are essential ingredi- 
ents for studies of relativistic heavy ion collisions at large 
net baryon densities. 

We find the observed non-monotonic behavior of T e ff 
with hadron mass to be a direct consequence of the non- 
equilibrium transport in the hadronic phase. Observa- 
tions like the centrality dependence of freeze out param- 
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eters and the difference in particle and anti-particle spec- 
tra and elliptic flow, on the other hand, can, at least par- 
tially, be explained by the conservation of baryon charge. 

We have argued that the observed difference in the 
elliptic flow between positively and negatively charged 
pions and kaons may likely be a consequence of the con- 
served net strangeness and isospin. Thus, future the- 
oretical studies of heavy ion collisions at low energies 
should take into account the explicit conservation of these 
charges in addition to the baryon number conservation. 
We further pointed out that a transport treatment of the 
final hadronic phase should include both proton - anti- 
proton annihilation and production processes to ensure 
detailed balance. This is essential to draw firm conclu- 
sions about the physics leading to subtle difference be- 
tween particle and anti-particle observables at low colli- 
sion energies, such as the observed difference in the ellip- 
tic flow. Finally, we have pointed out that local fluctua- 
tions of the baryon number may lead to a biased deter- 



mination of the event plane which may result in artificial 
differences between particle and anti-particle flow observ- 
ables. 
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